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Lithium-ion  battery  performance  is  affected  by  cut  surface  quality  during  the  electrodes’  cutting  process. 
Currently,  die  cutting  and  rotary  knife  slitting  have  been  used  to  cut  prismatic  and  cylindrical  electrodes, 
respectively.  Both  techniques,  which  require  expensive  tooling  that  wear  out  over  time,  result  in  process 
instability  and  poor  cut  quality.  These  will  cause  an  internal  short  circuit  and  significant  heat  generation 
in  the  cells.  Laser  cutting,  proved  and  widely  utilized  in  the  industry,  can  solve  the  abovementioned 
problems  by  improving  cut  surface  quality  due  to  many  advantages,  such  as  no  tool  wear,  high  energy 
concentration,  fast  processing  speed,  very  narrow  Heat  Affected  Zone,  applicability  to  nearly  all  mate¬ 
rials,  and  flexibility  of  laser  power.  Investigating  underlying  physical  phenomena  with  numerical  analysis 
provides  significant  advantages  to  fully  utilize  the  remote  laser  cutting  of  electrodes  for  lithium-ion 
batteries.  In  this  paper,  a  mathematical  model  of  three-dimensional  self-consistent  remote  laser 
cutting  is  presented  for  anode  (graphite-coated  copper)  of  lithium-ion  batteries.  Computational  and 
experimental  results,  obtained  by  using  laser  power  of  450  W  and  scanning  speed  of  5  m  s-1,  show  a  two- 
step  melt  pool  geometry  and  copper  composition  increase  on  the  material  interface. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

To  reduce  CO2  emissions,  the  United  States  Environmental 
Protection  Agency  (EPA)  set  up  stringent  emission  standards  in 
2009.  To  meet  these  stringent  standards,  automakers  have  focused 
heavily  on  the  development  of  Hybrid  Electrical  Vehicles  (HEVs), 
Electrical  vehicles  (EVs),  and  Plug-in  Hybrid  Electrical  Vehicles 
(PHEVs).  To  develop  these  vehicles  in  a  more  efficient  and  envi¬ 
ronmentally  friendly  manner,  an  understanding  of  automotive 
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battery  technology  is  essential.  Among  many  of  the  battery  tech¬ 
nologies,  the  lithium  ion  battery  represents  the  most  promising 
technology  due  to  its  high  energy-to-weight  ratio,  high  power 
density,  high  energy  density,  and  lighter  weight  than  other  energy- 
equivalent  secondary  batteries  [1,2]. 

Much  research  has  been  focused  on  improving  lithium  ion 
battery  performance:  developments  of  electrode  materials  and 
battery  manufacturing  processes  [3].  Many  researchers  have  been 
focused  on  development  of  electrode  materials.  Broussely  et  al.  [4] 
described  a  short  historical  perspective  of  positive  materials  by 
comparing  battery  performance  parameters  such  as  specific 
capacity.  Recently,  more  attention  has  been  given  to  cathode 
materials  to  improve  capacity  for  lithium  ion  batteries  through 
nanotechnology  [5-8].  Furthermore,  to  improve  the  capacity  of 
anode  materials,  silicon  or  silicon-based  materials  have  been 
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suggested  as  potential  alternatives  for  the  next  generation  of 
lithium  ion  battery  systems  [9-13]. 

Although  lithium  ion  battery  manufacturing  process  is  also  an 
important  issue,  little  research  has  been  done  in  this  area.  From  the 
manufacturing  process  of  prismatic  lithium  ion  cells,  several 
problems  have  been  reported  during  the  slitting  process.  Currently, 
die  cutting  and  rotary  knife  slitting  have  been  used  to  cut  the 
prismatic  and  cylindrical  electrodes,  respectively  [14-16].  Both 
techniques  require  relatively  expensive  tooling  that  wears  out  over 
time.  This  tool  wear  results  in  process  instability  and  requires 
additional  cost  and  downtime  to  replace  it.  Furthermore,  a  redesign 
of  mechanical  cutting  processes  demanding  extra  expense  is  an 
inevitable  step,  since  cell  specifications  vary  depending  on  the 
application,  such  as  electronic  devices  and  electric  vehicles.  The 
most  significant  disadvantage  due  to  the  tool  wear  is  the  problem  of 
inconsistent  cut  surfaces  with  defects,  which  are  characterized  by 
burrs,  delamination,  edge  bending,  and  micro-sized  material 
attachments,  as  shown  in  Fig.  1.  These  defects  cause  the  penetration 
of  the  separator,  giving  rise  to  internal  short  circuits.  Moreover,  the 
quality  of  a  poor  cut  surface  makes  electrodes  susceptible  to  higher 
electrical  stress,  which  can  lead  to  significant  heat  generation  and 
possible  thermal  runaway  [17].  Internal  short  circuits  and  signifi¬ 
cant  heat  generation  eventually  can  result  in  catastrophic  failure  of 
the  entire  system. 

Laser  cutting  is  proven  and  widely  utilized  in  industry  [18].  This 
highly  efficient  and  reliable  manufacturing  method  can  solve  the 
abovementioned  problems  since  it  has  many  advantages,  such  as 
no  tool  wear  due  to  contact-free  process,  high  energy  concentra¬ 
tion,  low  noise  level,  fast  processing  speed,  very  narrow  Heat 
Affected  Zone  (HAZ),  applicability  to  nearly  all  materials,  and 
flexibility  of  laser  power  [18].  Most  importantly,  laser  cutting  is 
relatively  easily  applicable  to  different  cell  specifications  with  less 


Fig.  1.  Micro-sized  material  attachments  (a),  edge  bending  (b),  and  inconsistent  cut 
quality  (c). 


cost,  since  the  redesign  of  laser  cutting  processes  is  relatively 
simpler  and  cheaper  compared  to  the  redesign  of  mechanical 
cutting  processes. 

While  laser  cutting  has  many  advantages,  the  application  of 
laser  cutting  to  size  the  lithium  ion  batteries  has  been  investigated 
by  few  researchers.  The  achievable  cutting  speeds  and  cut  qualities 
of  the  electrodes  for  the  lithium-ion  battery  have  been  studied  with 
a  pulsed  solid  state  laser  and  a  single  mode  fiber  laser  in  combi¬ 
nation  with  fixed  optics  and  2D  scanning  optics  [15].  High  speed 
laser  cuttings  of  electrodes  for  the  lithium-ion  battery  using  single 
mode  fiber  lasers  have  also  been  investigated  by  Patwa  et  al.  [16]. 
They  illustrated  the  achievable  highest  cutting  speed,  the  effect  of 
the  focus  beam  and  the  number  of  cutting  passes.  A  comparison 
between  a  Continuous  Wave  (CW)  laser  and  a  pulsed  laser  cutting 
of  electrodes  for  lithium-ion  batteries  has  additionally  been 
examined,  and  achievable  cutting  speeds,  qualities,  and  the  appli¬ 
cation  areas  for  the  two  operating  modes  have  been  also  discussed 
[19].  These  experimental  conclusions  have  been  derived  through 
a  trial  and  error  process  [19].  Investigating  underlying  physical 
phenomena  with  numerical  analysis  provides  significant  impacts  to 
fully  utilize  the  remote  laser  cutting  of  electrodes  for  lithium-ion 
batteries.  Since  the  understanding  of  physical  phenomena  for 
each  material  during  laser  cutting  is  an  essential  step  prior  to 
investigating  the  laser  cutting  of  electrodes,  numerical  studies  by 
Lee  and  Mazumder  [20]  of  laser  processing  parameters  on  current 
collectors  for  lithium-ion  batteries  have  been  done.  Thresholds  of 
laser  processing  parameters,  such  as  laser  power  and  scanning 
speed  for  the  laser  cutting  of  current  collectors,  have  been  found.  A 
mathematical  model  of  laser  cutting  of  an  anode  is  developed  in 
this  paper.  This  paper  is  organized  as  follows. 

First,  a  mathematical  model  of  three-dimensional  self-consis¬ 
tent  laser  cutting  is  depicted  for  anode  of  lithium-ion  batteries. 
Second,  an  experimental  set  up  and  procedures  are  stated.  Third, 
simulation  results  are  analyzed  in  terms  of  penetration  time,  depth, 
width,  and  absorptivity.  Furthermore,  melt  pool  flow,  melt  pool 
geometry,  and  copper  composition  changes  on  the  material  inter¬ 
face  are  discussed.  Fourth,  experimental  results  are  compared  with 
simulation  results.  Finally,  the  conclusions  drawn  from  this  work 
are  summarized. 

2.  Mathematical  modeling 

Laser  cutting  of  an  anode  material  includes  interconnected 
physical  phenomena.  When  a  laser  irradiates  a  substrate,  the 
substrate  absorbs  light  energy.  This  absorbed  energy  is  enough  to 
melt  and  evaporate  materials.  In  the  melted  area,  called  the  melt 
pool,  solid  and  liquid  materials  coexist  so  that  they  form  a  mushy 
zone.  Since  this  mushy  zone  consists  of  solid  and  liquid  phases, 
thermo-physical  material  properties  are  changed  based  on  the 
liquid  and  solid  fraction. 

In  the  evaporation  interface,  a  very  thin  layer  of  several  mean  free 
paths  is  observed.  This  thin  layer  is  called  a  kinetic  Knudsen  layer. 
Across  the  kinetic  Knudsen  layer,  steep  changes  of  temperature, 
pressure,  and  density  occur,  in  which  the  continuum  hypothesis 
fails.  Thus,  a  mathematical  discontinuity  needs  to  be  considered  as 
a  jump  condition  with  back  pressure  [21].  Furthermore,  mass  and 
thermal  energy  loss  on  the  interface  are  taken  into  account  in  this 
mathematical  model  through  the  jump  condition  [22]. 

In  addition  to  the  evaporation,  there  are  other  physical 
phenomena  that  can  be  observed  at  the  Liquid/Vapor  (L/V)  inter¬ 
face.  Vaporization  behaves  as  a  repulsive  force  on  the  melt  pool, 
and  is  referred  to  as  recoil  pressure.  This  recoil  pressure  generates 
a  narrow  and  deep  penetration  hole,  called  a  keyhole  or  deep 
penetration  hole.  This  keyhole  changes  the  geometry  of  the  L/V 
interface  and  absorbs  more  energy  by  internal  multiple  reflections 
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[23].  Moreover,  surface  tension  exists  due  to  the  geometry  change 
and  temperature  gradient. 

Different  from  single  materials,  electrodes  are  a  sandwiched 
composite  material.  Due  to  this  sandwiched  shape,  the  interface  of 
each  material  needs  to  be  carefully  considered.  A  material  interface 
of  the  two  sandwiched  materials  plays  a  significant  role  in  the  laser 
cutting  of  electrodes.  On  the  material  interface,  melted  materials 
are  mixed  and  compositions  change.  These  composition  changes 
vary  the  mixture  properties.  A  concentration  equation  is  used  to 
keep  track  of  composition  variation  on  the  boundary  of  the  current 
collector  and  active  electrode  materials  during  the  laser  cutting 
process.  Material  properties  of  the  mixture  on  the  boundary  of  the 
current  collector  and  active  electrode  materials  are  realized  with 
the  aid  of  the  phase  diagram  and  the  composition  equation. 

2.1.  Assumptions 

•  Laminar,  Newtonian,  incompressible  flow. 

•  Constant  material  properties. 

•  Graphite  has  no  liquid  phase. 

•  Laser  beam  has  Gaussian  distribution. 

•  Material  properties  between  vapor  and  liquid  are  smoothed 
out  for  numerical  purposes. 

•  Material  properties  for  high  temperature  range  are  extrapolated. 

•  Phase  diagram  of  copper  and  graphite  is  extrapolated. 


^2Ll  +  v-(puu)  =  -^  +  v-(mVu)  +  Su(x,y,z,t)  (2) 

-A- — \-v-(puv)  =  —  3^  +  V- (/iVv)  +  Sv(x,y,z,  t)  (3) 

+  V  •  (puw)  =  ~  +  V ■  (pWw)  +Sw(x,y,z,t )  (4) 

^(p^-pLCut)  /  —  \ 

— — ^ - —  +  u-v(pCpi_CuTJ  =  V-{kVT)+SE(x,y,z,t)  (5) 

^p  +  V-(u  pc)  =  V-(pDVc)  +Sc(x,y,z,  [)  (6) 

where  u  is  the  liquid  velocity  vector,  x  is  the  Cartesian  coordinate 
along  the  beam  scanning  direction,  y  is  the  width-wise  Cartesian 
coordinate,  z  the  depth-wise  Cartesian  coordinate,  p  is  the  density,  p  is 
the  viscosity,  p  is  the  pressure,  k  is  the  thermal  conductivity,  Cpl  Cu  is 
the  average-specific  heat  of  the  liquid,  and  t  is  the  time.  Sc,  Sm,  Su,  Sv,  Sw, 
SE,  and  Sc  are  the  source  terms  for  continuity,  x-momentum,  y- 
momentum,  z-momentum,  energy,  and  species  equations.  These 
source  terms  include  boundary  conditions  for  the  S/L  and  L/V  inter¬ 
faces.  These  source  terms  are  summarized  in  Table  1.  The  source  terms 
consist  of  the  boundary  conditions  of  the  L/V  and  S/V  interfaces. 


2.2.  Governing  equations 


2.3.  Source  terms  as  boundary  conditions  of  the  L/V  interface 


An  anode  material  consists  of  copper  as  a  current  collector  and 
graphite  as  an  active  electrode  material  as  shown  in  Fig.  2.  Since 
graphite  is  sublimated,  there  exist  two  phases,  solid  and  vapor,  at 
the  early  stage  of  the  simulation.  To  realize  this  sublimation  char¬ 
acteristic  of  graphite,  several  modifications  were  made.  Liquid 
viscosity  is  assumed  to  be  high  compared  to  the  gas  phase  since  the 
proposed  mathematical  model  solves  both  liquid  and  gas  phases. 
Moreover,  liquid  properties  are  set  to  be  the  same  as  solid  prop¬ 
erties,  and  surface  tension  is  ignored. 

However,  once  the  depth  of  a  deep  penetration  hole  reaches  the 
interface  of  graphite  and  copper,  the  liquid  phase  of  copper  and  the 
solid  phases  of  copper  and  graphite  exist  together.  When  this  mushy 
zone  is  formed,  composition  varies  at  the  interface.  Thus,  thermo¬ 
physical  properties  of  solid  and  liquid  mixtures  of  the  two  materials 
need  to  be  carefully  considered  based  on  the  composition  changes  and 
mass  fractions  of  the  two  materials.  To  incorporate  this  mushy  zone, 
this  study  adopted  a  method  developed  by  Bennon  and  Incropera  [24], 
which  is  an  extension  of  the  classical  mixture  theory.  Principles  of 
development  of  conservation  equations  are  as  follows.  First,  mixture 
components  may  be  viewed  as  isolated  subsystems  if  interactions  with 
other  mixture  components  are  properly  treated.  Second,  all  properties 
of  the  mixture  are  mathematical  consequences  of  the  component 
properties.  Finally,  the  mean  collective  mixture  behavior  is  governed 
by  equations  similar  to  those  governing  the  individual  components. 

V-u  =  Sm(x,y,z,t)  (1) 


Graphite  (90  pm) 
Copper  (10  pm) 


The  first  term  on  the  right  hand  side  of  the  source  terms  of 
continuity  equation  is  the  boundary  condition  taking  into  account 
evaporation  on  the  L/V  interface.  To  implement  vaporization,  jump 
conditions  with  back  pressure  are  used  to  take  into  account 
discontinuity  at  the  interface  caused  by  the  Knudsen  layer  [21  ].  The 
Knudsen  layer  implemented  as  jump  conditions  provides  the  net 
mass  flux,  rh"vap,  and  energy  flux  q"vap  [21,22,25].  The  second  terms 
on  the  right  hand  side  of  the  source  terms  of  momentum  equations, 
the  entire  term  within  parentheses,  are  boundary  conditions  for  the 
L/V  interface.  This  model  contains  both  the  capillary  and  thermo¬ 
capillary  forces,  similar  to  Ki  [26].  A  third  term  on  the  right  hand 
side  of  the  source  terms  of  the  energy  conservation  equation  is  the 
boundary  condition  of  the  L/V  interface. 

This  term  includes  energy  absorption  by  multiple  reflections, 
energy  losses  due  to  evaporation,  and  radiation  on  the  L/V  interface. 
Ql/v  is  the  spatial  laser  beam  distribution  after  multiple  reflections 
obtained  by  using  the  ray  tracing  method  since  total  energy 
transfer  from  the  laser  beam  to  a  work  piece  is  significantly 


Table  1 

Source  terms  for  mass,  x-momentum,  y-momentum,  z-momentum,  energy,  and 
species  equations. 

Equations  Source  term 


Mass 

Sm(x,y,z,  t) 

x-momentum 

Su{x,y,z,  t) 

y-momentum 

Sv(x,y,z,  t) 

z-momentum 

Sw(x,y,z,  t) 

Energy 

SE(x,y,z,  t) 

Species 

Sc(x,y,z,  t) 

= 

=  "Y  ^u-ex'^n(V-n)-Vsr^5(0) 
=  -^^-ey^n(V-n)-VsT^5(0) 

=  “f^w-ez'(ffn(v'n)"Vs  Twf)5^ 

,  d(pfi)  ,  a(p/sACpr) 
at  at 

+  [tf/V  -  Qevap  -^(T4-7l)]5((/») 

=  V-  [pDV(q  —  c)]  —  V  -  [p(C\  —  c)u] 


Fig.  2.  Anode  material. 
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increased  due  to  the  multiple  reflections  [23].  q"vap  is  energy  loss 
due  to  evaporation  and  ae(T4  -  T4  )  is  energy  loss  due  to  radiation, 
cr  and  e  are  Stefan-Boltzmann  constant  and  emissivity,  respectively. 
Note  that  the  simultaneous  calculation  of  the  vapor  phase  and  the 
liquid  phase  allows  for  the  exclusion  of  convection  heat  loss  [26]. 

The  formula  for  the  multiple  reflections  is  implemented  with 
ray  tracing  method  [27].  The  material  properties  are  smoothed 
from  the  liquid  to  the  vapor  phase  since  the  L/V  interface  has 
significantly  different  properties  [27].  Using  ray  tracing  method,  the 
laser  energy  is  absorbed  into  the  wall  at  each  reflection.  At  each 
reflection,  the  absorption  coefficient  on  the  L/V  interface  is  ob¬ 
tained  as  a  summation  of  the  two  absorption  coefficients,  consid¬ 
ering  the  mass  fractions,  and  can  be  written  as 

A{a)  =  (1 -/cs)A0_cu(cosa)‘J+/c/o_c(cosa;)<3  (7) 


h  =  fcuMcu,  +/cu|hcu,  +/c7!c  (14) 

The  relationship  between  the  volume  and  mass  fractions  is 


f  PCus§Cus 

/Cus  =  p 

(15) 

r  PCu\gCu\ 

/cui  _  P 

(16) 

c  PcsgCs 

fc>=  p 

(17) 

The  mass  and  volume  fractions  must  add  to  unity  as  follows 
fcus  +/cs  +/cu!  =  1  (18) 


where  a  is  the  incident  angle,  /Cs  is  solid  mass  fraction  of  graphite, 
A)_cu  is  the  absorption  coefficient  of  copper  on  the  flat  surface,  A0  c 
is  the  absorption  coefficient  of  graphite  on  the  flat  surface,  and  q  is 
the  characteristic  of  the  material  [28].  The  calculation  of  the 
multiple  reflections  is  terminated  when  the  rays  escape  the 
computational  domain. 


2.4.  Source  terms  as  boundary  conditions  of  the  S/L  interfaces 

The  boundary  conditions  of  the  S/L  interface  are  carefully 
considered  based  on  the  composition  changes  and  mass  fractions  of 
the  two  materials.  To  incorporate  this,  this  study  adopted  a  method 
developed  by  Bennon  and  Incropera  [24].  The  first  terms  on  the  right 
hand  side  of  the  source  terms  of  momentum  equations  are  the 
darcian  damping  force  to  treat  the  mushy  region  as  a  fine  permeable 
solid  matrix.  /i\  is  the  liquid  viscosity,  and  p\  is  the  liquid  density.  I<  is 
the  isotropic  permeability  and  is  assumed  to  vary  with  liquid  volume 
fraction  according  to  the  Kozeny-Carman  equations  [29]: 


I< 


I<r 


S\ 


(!-&)■ 


(8) 


The  first  and  second  terms  on  the  right  hand  side  of  the  source 
terms  of  the  energy  and  species  equations  are  the  expression  of  the 
mushy  zone.  These  are  achieved  from  Bennon’s  continuum  model 
[24]  by  linearizing  the  phase  enthalpies.  Derivations  are  briefly 
summarized  here.  Defining  the  mass  fraction,  f  and  the  volume 
fraction,  g,  with  phases  and  materials,  the  velocity,  density,  thermal 
conductivity,  mass  diffusion  coefficient,  concentration,  and 
enthalpy  for  the  liquid  and  solid  mixtures  of  the  two  materials  are 

U  =  /cusuCUs  +/culuCU|  +/csuCs  (9) 


P  —  gcusPcus  +  gcu1Pcul  +£csPcs 

/<  =  (pH,+pL+scL y1 

V^CUs  fccil!  kcs) 


(10) 

(11) 

(12) 


(13) 


&cus  +  goi!  +  &cs  =  1  (19) 

and  the  solid  and  liquid  mass  fractions  are 


fs  =  fcus  +/cs  (20) 

fi  =  fcui  (21) 

These  defined  mixture  variables  are  considered  in  the  continuity 
equation,  momentum  equation  in  x-,  y-,  and  z-  directions,  energy 
equations  and  species  equations.  The  phase  enthalpies  are  obtained 
as 


T 

hCUs  =  J  Cps_CudT  (22) 

0 

T 

K  =  J  Cps_cdT  (23) 

0 

7m  T 

l^CU]  —  J  Cps.cudr  +  Lm  +  J  Cpi_CudT  (24) 

0  Tm 

where  Tm  is  the  melting  temperature,  Lm  is  the  latent  heat  of  fusion, 
CpLcu  is  the  constant-pressure  specific  heat  of  the  liquid  phase  of 
copper,  Cps  cu  is  the  constant-pressure  specific  heat  of  the  solid 
phase  of  copper,  and  Cps_c  is  the  constant-pressure  specific  heat  of 
the  solid  phase  of  graphite.  Furthermore,  the  average  specific  heats 
of  solid  and  liquid  (Cps  and  Cpl)  are  defined  as 

T 

Cps_Cu  —  j  J  Cps_Cudl'  (25) 

o 


(-pl_Cu  —  j  _j  J  ^pl-Cud^ 


(26) 


6ps_c  —  j  I  6ps_cdT 


D  =  fcuP 

C  =  fcu^Cui  +  /cusQlus 


0 


T 


(27) 
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=  Cps_cT  (28) 

^Cus  —  Cps_Cu^  (29) 

^CU[  —  Cps  Cu^m  +  +  CpLCuC^  ~  Fm) 

—  CpLCu^  +  (Cps_Cu  —  Cpicu^m  +  =  CpLCu^  +  L 

(30) 

where  L  =  (Cps_Cu  -  Cpl  Cu)Tm  +  Lm.  By  using  the  defined  mixture 
variables,  the  boundary  conditions  for  the  S/L  interface  can  be 
calculated  by  substituting  Eqs.  (28)— (30)  into  Eq.  (14). 

^  —  (fcsCpsJ2  3" /cus^-ps_Cu  H~ /cui^pLCu^ ^  ~F  fc\i\L  =  CpT  +  fcUlL 

(31) 

where  Cp  is  the  average  mixture  specific  heat,  defined  as 

(-p  =  /cs  (-ps_C  +  fc  us  ^-ps_Cu  +  /cui  ^pl_Cu  (32) 

and  can  be  written  as 

Cp  =  Cpi  qu  — /sACp  (33) 

ACp  =  ^  (cpLCu  -  Cps  Cu)  +  ^  (Cpi.cu  -  Cps_c)  (34) 

Now,  the  first  and  second  terms  on  the  right  hand  side  of  the 
source  terms  of  the  energy  equation  can  be  obtained  by 
substituting  the  above  relations,  as  described  by  Ki  [30].  At  this 
point,  the  mass  fractions  have  to  be  obtained  and  updated  for  each 
numerical  iteration.  The  phase  diagram  can  provide  these  values. 
Since  the  currently  available  phase  diagram  provides  limited 
information,  the  rest  of  the  phase  diagram  is  extrapolated  to  the 
sublimation  temperature  of  pure  graphite  using  piecewise  contin¬ 
uous  functions  as  shown  in  Fig.  3.  The  solid  and  liquid  fractions  can 
be  obtained  by  the  lever  rule  (Fig.  3). 

A  boundary  condition  for  the  S/L  interface  of  the  composition 
equation  can  be  written  as 

=  1  (35) 

where  c0  is  initial  composition  of  copper  and  f  indicates  the  S/L 
interface.  Considering  the  boundary  condition  during  the  compu¬ 
tation  on  the  S/L  interface,  the  concentration  of  copper  into 
graphite  can  be  fully  included  in  this  mathematical  model. 


Fig.  3.  Extended  C-Cu  phase  diagram. 


2.5.  Level  set  equation  for  the  L/V  interface  tracking 

During  laser  cutting,  the  L/V  interface  is  continuously  changing 
due  to  complex  physical  phenomena  at  the  interface  such  as 
evaporation,  convective  flow,  multiple  reflections,  energy  absorp¬ 
tion,  energy  loss,  and  significantly  different  properties  between 
liquid  and  vapor.  This  section  describes  the  implementation  of 
these  interfacial  phenomena  into  the  proposed  model.  A  key  factor 
in  the  investigation  of  the  interface  shape  and  the  process  physics  is 
the  tracking  of  the  free  surface.  Osher  and  Sethian  [31]  developed 
the  level  set  method.  This  method  has  increasingly  been  used  for 
many  complex  problems  due  to  its  easy  implementation,  straight¬ 
forward  concept,  ability  to  handle  surface  merging  and  separation 
automatically,  and  easy  achievements  of  geometric  qualities,  such 
as  the  surface  normal  and  the  curvature  on  the  surface. 

The  L/V  interface,  which  is  the  surface  of  interest,  is  set  as  the 
zero  iso-surface  or  </>(x,t),  to  implement  the  level  set  method.  By 
equating  this  zero  iso-surface  with  zero  and  taking  the  material 
derivative,  we  can  get  the  level  set  equation  in  terms  of  a  hyperbolic 
partial  differential  equation. 

^  +  FW|=0  (36) 

where  F  is  the  speed  function.  The  detail  derivation  of  this  equation 
can  be  found  in  Ref.  [32]. 

During  the  laser  cutting  process,  the  L/V  interface  in  the  level  set 
equation,  i.e.  </>  =  0,  is  evolved  on  the  basis  of  the  speed  function, 
which  is  composed  of  vaporization  and  fluid  flow.  By  solving  the 
proposed  momentum  equation,  the  fluid  flow  can  be  obtained.  To 
implement  vaporization  into  the  speed  function,  jump  conditions 
with  back  pressure  are  used  to  take  into  account  a  thin  disconti¬ 
nuity  of  the  order  of  a  molecular  mean  free  path  thickness  at  the 
interface.  Hence,  a  final  form  of  the  speed  function  is 

F  —  ^evap  +U  n  (37) 

where  u  is  the  liquid  velocity  vector,  n  is  the  normal  vector,  and 
Fevap  is  the  vaporizing  mass  flux,  which  is  assumed  to  be  perpen¬ 
dicular  to  the  L/V  interface  and  obtained  by  the  Mott-Smith-type 
solution  [21,22,25]. 

3.  Experiments 

Experimental  investigation  is  designed  with  an  IPG  single  mode 
CW  fiber  laser  (model:  YLR-500SM)  with  the  diameter  of  11  pm  and 
the  wavelength  of  1.070  pm.  Its  maximum  laser  power  is  500  W.  The 
laser  beam  is  delivered  by  a  fiber  with  a  10  pm  core-fiber  diameter 
and  collimated  to  a  beam  diameter  of  13  mm.  This  fiber-delivered 
laser  beam  is  deflected  and  moved  by  a  2D  galvo-scanner,  made 
by  Scanlab  (model:  SK1020K),  at  high  dynamic  speeds  with 
maximum  speeds  up  to  5  m  s-1.  The  2D  galvo-scanner  is  connected 
to  an  80  mm  F-theta  lens,  which  keeps  the  scanning  image  field  flat 
within  a  35  x  35  mm  square  field  size.  The  laser  beam  is  measured 
using  FocusMonitor,  manufactured  by  Primes.  The  measured  spot  is 
symmetrical  and  has  closely  Gaussian  energy  distribution  as  shown 
in  Fig.  4.  At  the  focus  position  the  measured  value  of  spot  size  is 
approximately  11  pm.  The  Rayleigh  length  of  70  pm  and  M2  value 
(86%)  of  1.3  are  obtained  from  this  measurement. 

Fig.  5  shows  the  schematic  of  the  experimental  set  up.  Cutting  is 
performed  as  the  laser  beam  is  moving  across  a  stationary  anode 
surface  using  a  high  speed  galvo-scanner.  A  vacuum  fixture  is  used 
to  hold  anodes  and  prevent  any  movement  during  the  entire 
cutting  process.  A  narrow  groove,  which  is  less  than  0.5  mm,  was 
placed  in  the  top  plate  of  the  fixture  where  the  cutting  takes  place. 
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Fig.  4.  Measured  laser  beam  shape. 


This  prevents  any  attachments  of  melt  that  moves  out  of  cut  kerf  to 
the  fixture  plate.  The  fixture  is  mounted  on  the  horizontal  plane  of 
a  CNC  motion  system  and  the  2D  galvo-scanner  is  mounted  on  its 
vertical  axis.  Using  this  precise  3  axis  positioning  system,  the 
scanner  was  aligned  with  the  fixture  and  placed  at  the  correct 
vertical  distance  from  the  anode.  There  is  no  process  gas  applied 
during  cutting.  Cutting  is  performed  in  a  clean  and  dust  free 
environment  with  suitable  fume  exhaust  systems.  The  single  side- 
coated  anode  is  used  as  seen  in  Fig.  2.  The  anode  is  placed  on  the 
fixture  as  graphite  faces  upward  and  the  laser  beam  irradiates  on 
the  graphite  surface.  Thickness  of  graphite  and  copper  are  90  pm 
and  10  pm,  respectively.  Average  article  size  of  graphite  is  20.36  pm. 
The  density  of  the  graphite  part  of  the  coated  anode  is  120  g  m-2. 
Graphite  is  composed  of  active  graphite  powder  (92.5%),  CMC 
(1.5%),  SBR  liquid  (2.5%),  and  Conductive  (3.5%). 

4.  Results  and  discussions 

4.2.  Simulation  preparations  and  solution  schemes 

A  CW  Gaussian  laser  beam  is  used  and  scanned  in  the  positive  x 
direction.  A  focused  beam  diameter  is  11  pm.  It  is  assumed  that  the 
beam  profile  has  no  z  dependency. 


2D  gavlo-scanner 


Vacuum  hose 


/(r)  =  /0exp  - 


2  r2 


(38) 


where  fib  is  the  beam  radius  and  Io  is  the  centerline  laser  intensity. 
The  laser  power  and  scanning  speed  chosen  for  this  study  are  450  W 
and  5  m  s*1,  respectively.  The  computational  domain  of  the  simula¬ 
tion  is  55  pm  x  30  pm  x  130  pm.  Non-uniform  mesh  is  used  to  save 
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Fig.  5.  Experimental  set  up. 


Fig.  6.  Simulation  flow  diagram. 
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Table  2 

Material  properties  of  copper. 


Property 

Value 

Melting  temperature 

1357.77  (K) 

Normal  boiling  temperature 

2835.15  (K) 

Critical  point  temperature 

8280  (K) 

Liquid  density 

7920  (kg  mr3) 

Solid  density 

8960  (kg  mr3) 

Kinematic  viscosity  [36] 

3.50E  -  07  (m2  s’1) 

Surface  tension  [37] 

1.257  -  0.0002*(T  -  1356) 

(N  ITT1) 

Latent  heat  of  vaporization 

5.23E  +  06  (J  kg"1) 

Latent  heat  of  fusion 

2.05E  +  05  (J  kg"1) 

Solid  thermal  conductivity  [38] 

317  (W  mr1  K"1) 

Liquid  thermal  conductivity  [38] 

157  (W  m"1  K"1) 

Liquid  constant-pressure  specific  heat 

571.6218  (J  kg"1  K"1) 

Solid  constant-pressure  specific  heat  [39] 

385  y  kg-1  K-1) 

Liquid  thermal  diffusivity 

3.62E  -  05  (m2  s"1) 

Solid  thermal  diffusivity 

7.63E  -  05  (m2  s"1) 

Laser  absorptivity  for  flat  surface 

0.05 

computational  time.  The  smallest  discrete  mesh  distance  is  chosen  as 
0.6  pm.  A  staggered  grid  is  utilized  to  attain  the  physically  proper 
pressure  field  [33].  One  side  graphite-coated  copper  is  chosen  as  the 
substrate  material,  where  the  thickness  of  the  anode  is  100  pm. 
Graphite  of  90  pm  and  copper  of  10  pm  are  placed  as  the  top  and 
bottom  layers  of  the  anode,  respectively.  The  properties  of  these 
materials  are  shown  in  Tables  2  and  3.  An  implicit  finite  difference 
method  is  used  for  the  discretization  of  the  governing  equations.  To 
solve  the  matrix  equation  and  coupled  pressure-velocity  fields  on 
the  staggered  grid  system,  the  conjugated  Gradient  Stabilized 
(CGSTAB)  method  and  Semi-Implicit  Method  for  Pressure-Linked 
Equation-Consistent  (SIMPLEC)  are  adopted,  respectively.  In  addi¬ 
tion,  a  second-order  space  convex  scheme  [34]  is  used  to  discretize 
the  spatial  derivative  of  the  level  set  equation.  Finally,  the  value  of  c 
for  all  of  the  material  properties  is  taken  0.8  and  7  x  min(  Az)  is  chosen 
for  £.  Az  is  the  shortest  length  of  the  smallest  mesh. 

A  flow  diagram  of  the  solution  procedure  in  one  time  step  is 
shown  in  Fig.  6.  First,  simulation  parameters  and  material  properties 
are  assigned.  In  addition,  the  jump  conditions  for  the  Knudsen  layer 
are  tabulated  as  an  initialization  step.  Second,  the  multiple  reflections 
are  executed.  Third,  the  temperature,  velocities,  and  concentration  of 
copper  are  solved  sequentially  in  an  inner  loop  with  coupled  solu¬ 
tions  until  the  coupled  solutions  satisfy  computational  residuals.  A 
maximum  iteration  is  set  to  prevent  unnecessary  computations  due 
to  slow  convergence,  even  if  the  residual  values  are  not  satisfied. 
Fourth,  the  level  set  equation  is  calculated  with  the  obtained  values  of 
evaporation  and  convection,  which  are  main  sources  of  interface 
deformation.  Finally,  the  mixture  properties  are  updated  with  a  given 
liquid  mass  fraction  through  the  phase  diagram.  This  procedure  is 
repeated  until  it  satisfies  the  designated  computation  time. 

4.2.  Characteristics  of  the  anode  laser  cutting 

In  the  following  sections,  the  top  layer  (90  pm)  is  called  the 
graphite  region;  the  bottom  layer  ( 10  pm)  is  called  the  copper  region. 


Table  3 

Material  properties  of  graphite. 


Property 

Value 

Sublimation  temperature  [40] 

4800  (K) 

Critical  point  temperature  [41  ] 

7811  (K) 

Solid  density  [42] 

1730  (kg  m"3) 

Latent  heat  of  sublimation  [42] 

5.98E  +  07  (J  kg"1) 

Solid  thermal  conductivity  [42] 

18.1  (W  m"1  K"1) 

Solid  constant-pressure  specific  heat  [42] 

2092.48  (J  kg  1  K"1) 

Solid  thermal  diffusivity  [42] 

5.00E  -  06  (m2  s"1) 

Laser  absorptivity  for  flat  surface 

0.81 

Penetration  depth,  width,  and  absorptivity  with  evolution  of  time  are 
shown  in  Fig.  7.  A  deep  penetration  hole  reaches  a  material  interface 
between  graphite  and  copper  at  0.1646  ps.  Around  this  material 
interface,  several  interesting  overall  behaviors  are  observed.  The 
slope  of  the  depth  increase  of  the  deep  penetration  hole  changes  over 
time.  Before  the  deep  penetration  hole  reaches  the  material  interface, 
its  slope  of  the  depth  increase  is  722.16  m  s-1.  Once  the  depth  of  the 
deep  penetration  hole  reaches  the  material  interface,  the  slope 
decreases  to  288.00  m  s-1.  The  width  increases  gradually  and  it 
becomes  31.196  pm  when  the  deep  penetration  hole  reaches  the 
material  interface.  This  gradual  increase  of  width  is  possible  since  the 
width  is  measured  on  the  top  surface  of  the  anode.  At  the  beginning  of 
the  simulation,  the  absorptivity  is  0.81,  which  is  the  absorption 
coefficient  of  graphite  as  a  flat  substrate.  The  maximum  absorptivity 
is  0.9588  at  which  point  the  deep  penetration  hole  reaches  the 
material  interface.  After  obtaining  the  maximum  value  at  0.1646  ps, 
absorptivity  starts  to  decrease  gradually  with  fluctuation  and  then  it 
increases  again  before  it  has  completed  thorough  cutting. 

Characteristics  of  penetration  depth  and  absorptivity  change 
significantly  at  the  point  where  the  deep  penetration  hole  reaches 
the  material  interface.  These  observations  are  closely  related  to 
each  other  and  due  mainly  to  the  composition  change  at  the 
material  interface.  An  absorption  coefficient  changes  its  value 
depending  on  the  mass  fraction  on  the  surface,  since  the  absorption 
coefficient  is  obtained  by  a  summation  of  the  two  absorption 
coefficients  of  the  graphite  and  copper  with  mass  fractions  as 
described  in  Eq.  (8).  Increasing  the  copper  concentration  around 
the  material  interface  contributes  to  decreasing  the  absorption 
coefficient  of  the  mushy  zone  so  that  absorptivity  decreases.  The 
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Fig.  7.  Penetration  depth  (a),  width  (b),  and  absorptivity  (c)  with  time. 
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Fig.  8.  Temperature  distribution  and  melt  pool  geometry  of  anode  (laser  power:  450  W,  scanning  speed:  5000  mm  s  at  (a)  0.043  ps,  (b)  0.066  ps,  (c)  0.126  ps,  (d)  0.164  ps,  (e) 
0.172  ps,  (f)  0.183  ps,  (g)  0.192  ps,  (h)  0.200  ps. 


decrease  of  absorptivity  provides  less  energy  absorption  on  the  L/V 
interface.  Therefore,  the  rate  of  depth  increases  changes. 

4.3.  Temperature  distribution  and  melt  pool  geometry 

Temperature  distribution  and  melt  pool  geometry  are  shown  in 
Fig.  8.  At  the  beginning  of  the  simulation,  the  deep  penetration  hole 
shows  a  very  smooth  and  clean  surface  as  seen  in  Fig.  8  (a)-(c).  The 
surface  shows  clearly  the  characteristic  of  the  laser  cutting  of 
graphite,  which  was  investigated  in  the  study  of  laser  cutting  on  an 
active  electrode  material  [35].  No  liquid  phase  exists  for  graphite 
since  it  sublimates.  Hence,  the  resulting  surface  is  smooth  and 
clean.  In  addition,  high  temperatures  are  observed  in  the  front  and 
back  side  of  the  L/V  interface  due  to  multiple  reflections.  When  the 
deep  penetration  hole  reaches  the  material  interface,  a  melt  pool 
starts  to  form  as  seen  in  Fig.  8  (d).  After  this  initial  melt  pool 
formation,  the  melt  pool  develops  and  shows  the  two-level  surface, 
two-step-like  shape,  as  shown  in  Fig.  8  (g).  This  two-level  surface 
can  be  explained  by  the  different  material  properties  between 


graphite  and  copper.  While  no  surface  changes  in  the  graphite 
region  due  to  its  sublimation  property,  melt  pool  flow  observed  in 
the  copper  region  changes  the  geometry  of  the  bottom  part  of  the 
deep  penetration  hole.  Therefore,  two-level  surfaces  form. 

Penetrations  on  the  side  and  front  of  the  deep  penetration  hole 
at  the  moment  of  full  penetration  are  observed  as  shown  in  Fig.  8 
(h).  The  movement  of  the  liquid  copper  causes  the  fluctuation 
motion  of  the  deep  penetration  hole  and  creates  an  uneven  melt 
pool  shape  in  the  copper  region.  This  uneven  shape  reflects  a  laser 
beam  toward  the  graphite  region  and  graphite  sublimates  due  to 
the  absorption  of  reflected  laser  energy.  Hence,  the  reflected  laser 
beam  forms  the  front  and  side  penetrations  and  changes  the  melt 
pool  geometry  dramatically  around  the  material  interface. 

4.4.  Copper  composition  distribution  and  melt  pool  flow 

Figs.  9  and  10  show  copper  composition  distribution  and  melt 
pool  flow,  respectively.  No  copper  composition  and  melt  pool  flow 
are  shown  in  Figs.  9  (a)-(c)  and  10  (a)-(c)  since  the  deep 
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Fig.  9.  Copper  composition  distribution  (laser  power:  450  W,  scanning  speed:  5000  mm  s  at  (a)  0.043  [is,  (b)  0.066  [is,  (c)  0.126  ps,  (d)  0.164  ps,  (e)  0.172  ps,  (f)  0.183  ps,  (g) 
0.192  ps,  (h)  0.200  ps. 


penetration  hole  forms  in  the  graphite  region.  When  the  deep 
penetration  hole  reaches  the  material  interface,  energy  transfers 
from  the  graphite  surface  to  the  solid  copper.  The  copper  increases 
its  temperature  to  the  melting  point  so  that  the  copper  starts  to 
melt  and  the  melt  pool  forms.  Next,  the  copper  concentration  starts 
to  increase  at  the  material  interface,  as  seen  in  Figs.  9  (d)  and  10  (d). 
At  the  beginning  of  the  mass  transfer,  higher  values  of  the  copper 
concentration  are  observed  on  the  front  and  side  walls  of  the  deep 
penetration  hole  as  shown  in  Fig.  9  (e).  In  Fig.  9  (f)-(h),  as  the  deep 
penetration  hole  is  becoming  deeper,  the  copper  concentration  is 
distributed  more  uniformly  on  the  copper  region;  copper  concen¬ 
tration  on  the  graphite  region  increases  around  the  deformed 
graphite  surface.  However,  the  back  side  of  the  deep  penetration 
hole  in  the  copper  region  shows  less  copper  concentration.  The 
front  melt  pool  provides  stronger  melt  pool  flow,  as  seen  in  Fig.  10 
(e)-(h).  Furthermore,  melt  pool  flow  is  observed  on  the  graphite 
region  due  to  the  existence  of  the  liquid  copper.  When  a  full 
penetration  occurs,  as  shown  in  Fig.  9  (h),  the  copper  concentration 
is  relatively  low  around  the  bottom  of  the  penetrated  hole.  Stronger 
flow  is  observed  around  the  deformed  graphite  surface,  as  shown  in 
Fig.  10  (h). 


4.5.  Comparison  between  simulation  and  experiments 

The  top  view  of  a  laser  cut  anode  is  shown  in  Fig.  11  (a).  A 
laser  irradiates  the  graphite  surface  as  shown  in  Fig.  2.  Kerf 
widths  on  the  top  and  bottom  layers  are  97.60  pm  and  27.90  pm, 
respectively.  Luetke  et  al.  [19]  defined  the  clearance  width  of  the 
metal  foil  of  anodes  for  lithium-ion  batteries,  which  is  the 
ablation  width  of  the  upper  coating.  This  study  also  observes  the 
clearance  width  of  29.10  pm.  These  experimental  results  clearly 
show  the  two-level  surface  on  the  material  interface,  which  is 
also  observed  in  the  simulation  result.  The  kerf  width  obtained 
by  simulation  is  31.19  pm  and  shown  in  Fig.  11  with  a  dashed 
line.  The  simulation  result  overestimates  the  kerf  width  of  the 
bottom  layer  by  7.18%,  which  is  in  good  agreement  with  the 
experimental  measurement.  However,  the  kerf  width  of  the  top 
layer  shows  discrepancy  compared  to  the  kerf  width  obtained  by 
the  simulation.  The  discrepancy  of  the  kerf  width,  which  is 
characterized  by  the  clearance  width,  can  be  explained  from  the 
computational  analysis.  The  fluctuation  motion  of  the  deep 
penetration  hole  due  to  the  liquid  flow  of  copper  creates  uneven 
melt  pool  shape  at  the  tip  of  the  deep  penetration  hole.  The 
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Fig.  10.  Melt  pool  flow  of  anode  (laser  power:  450  W,  scanning  speed:  5000  mm  s  ^  at  (a)  0.043  \is,  (b)  0.066  ps,  (c)  0.126  ps,  (d)  0.164  ps,  (e)  0.172  ps,  (f)  0.183  ps,  (g)  0.192  ps,  (h) 
0.200  ps. 


reflected  laser  beam  toward  the  graphite  region  from  the  uneven 
melt  pool  shape  sublimates  graphite.  Hence,  the  kerf  widths  of 
the  graphite  region  and  copper  region  are  different.  The  side 
view  of  the  laser  cut  anode  is  shown  in  Fig.  11  (b).  It  shows  good 
quality  of  cut  surface.  No  delamination,  edge  banding,  burrs,  or 
micro-size  attachments  on  the  material  interface  are  observed. 
At  the  bottom  of  the  anode,  debris  caused  by  re-solidification  is 
observed.  Copper  concentrations,  which  are  obtained  by  both 
the  simulation  and  experiment,  are  analyzed  along  the  perpen¬ 
dicular  line  as  shown  in  Fig.  12  and  its  values  are  plotted  in  Fig. 
13.  The  copper  concentration  of  the  experimentally  obtained 
anode  sample  was  obtained  through  Energy-Dispersive  X-ray 
spectroscopy  (EDX)  analysis. 

From  the  top  surface  to  the  depth  of  71.31  pm,  the  experi¬ 
mentally  measured  copper  concentration  is  lower  than  10%  and 
shows  no  trend.  The  copper  concentration  increases  significantly 
from  the  depth  of  71.31  pm  to  82.55  pm  and  reaches  its  value  of 
99.68%  at  the  depth  of  82.55  pm.  As  it  become  deeper,  the  copper 
concentration  gradually  decreases  from  the  99.68%  to  16.12%,  from 
the  depth  of  82.55  pm  to  the  material  interface,  and  then  it 
increases  again  and  becomes  100%  copper  concentration  at  the 


depth  of  100  pm.  For  the  computationally  obtained  copper 
concentration,  there  is  no  copper  concentration  observed  from  the 
top  surface  to  the  depth  of  73.96  pm.  From  the  depth  of  73.96  pm 
to  74.56  pm,  a  small  amount  of  copper  concentration  is  found.  The 
copper  concentration  increases  significantly  and  becomes  100%  at 
the  depth  of  84.73  pm.  The  copper  concentration  stays  at  100% 
until  at  the  depth  of  89.52  pm  and  then  decreases  up  to  78.02%  at 
the  depth  of  99.10  pm.  Both  experimentally  and  computationally 
obtained  copper  concentrations  share  three  characteristics.  First, 
the  copper  concentration  starts  to  increase  its  value  at  least  10  pm 
above  the  material  interface.  Second,  it  increases  very  sharply  up 
to  100%  copper  concentration  after  its  initial  increase.  Third,  the 
copper  concentration  decreases  once  100%  copper  concentration  is 
achieved.  The  reason  of  the  99-100%  copper  concentration  at 
about  10  pm  above  the  material  interface  is  the  direction  of  liquid 
copper  as  we  see  in  Fig.  10.  The  thorough  penetration  may  affect 
the  decrease  of  the  copper  concentration  near  the  material  inter¬ 
face.  In  addition,  the  flow  direction  removes  liquid  copper  from 
the  bottom  of  the  deep  penetration  hole  to  the  outward  direction. 
This  also  explains  the  99-100%  copper  concentration  at  about 
10  pm  above  the  material  interface. 


378 


D.  Lee  et  al.  /  Journal  of  Power  Sources  240  (2013)  368—380 


Fig.  11.  Top  view  (a)  and  side  view  (b)  of  laser  cut  anode. 


Although  both  computational  and  experimental  results  share 
similar  trends,  there  is  a  discrepancy  in  terms  of  kerf  widths  and 
the  copper  concentration.  The  discrepancy  may  be  caused  by 
three  main  reasons.  First,  graphite  is  mixed  with  additives  and 
binders  during  the  graphite  coating  process.  Hence,  its  material 
properties  could  be  different  from  the  material  properties  of  pure 
graphite,  adapted  in  this  simulation.  Second,  the  extended  copper 
and  graphite  binary  diagram  may  give  rise  to  this  discrepancy. 
The  liquid  fraction  can  be  overestimated  since  the  binary  diagram 
is  extended.  Finally,  the  simulation  does  not  take  into  account  the 


Copper  concentration/  % 


Fig.  13.  Comparison  of  the  copper  composition  between  experiment  and  simulation. 


reaction  with  oxygen,  or  the  air  suction  to  hold  the  sample 
surface  flat. 

As  mentioned  above,  the  properties  of  the  coated  graphite  layer 
are  different  from  those  of  pure  graphite.  The  graphite  layer  of 
experimental  electrode  probably  has  a  lower  solid  density,  a  higher 
porosity  and  a  lower  thermal  conductivity.  These  differences  may 
affect  entire  simulation  results.  First  of  all,  the  kerf  width  and  the 
slope  of  depth  penetration  may  be  decreased  since  the  coated 
graphite  layer  has  a  lower  solid  density,  a  higher  porosity  and 
a  lower  thermal  conductivity.  Furthermore,  flow  speed  and  pattern 
of  the  liquid  copper  could  be  changed  since  the  geometry  of  the 
first  layer  highly  influence  the  transport  phenomena  of  the  liquid 
copper  during  the  laser  material  interaction. 

The  simulation  analysis  in  this  study  found  strong  outward 
liquid  flow  at  the  front  of  the  deep  penetration  and  composition 
changes  on  the  graphite  region.  The  strong  outward  flow  of  liquid 
copper  gives  rise  to  composition  change  and  sharp  edge  formation 
on  the  cut  surface.  The  composition  change  on  the  material  inter¬ 
face  may  result  in  non-uniform  electrochemical  reactions  and 
higher  electrical  stress,  which  can  lead  to  significant  heat  genera¬ 
tion  and  possible  thermal  runaway  during  the  use  of  the  lithium- 
ion  battery  cell.  In  addition,  the  sharp  edge  could  penetrate  the 
separator  so  that  an  internal  short  circuit  occurs. 


Fig.  12.  Copper  composition  variation  along  the  line  both  for  experimental  sample  (a)  and  simulation  (b). 
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5.  Conclusions 


Nomenclature 


Laser  cutting  allows  an  improved  quality  of  cut  surface  and 
cutting  speed  during  lithium-ion  battery  manufacturing  processes. 
The  advantages  of  laser  cutting  can  be  maximized  by  under¬ 
standing  the  underlying  physics  during  the  laser  cutting  of  elec¬ 
trodes  for  lithium-ion  batteries.  To  understand  the  underlying 
physics,  a  mathematical  model  of  three-dimensional  self-consis¬ 
tent  laser  cutting  of  an  anode  for  lithium-ion  batteries  is  developed 
including  multiple  physical  phenomena  such  as  heat  transfer, 
composition  changes,  fluid  flow,  recoil  pressure,  capillary  and 
thermo-capillary  forces,  multiple  reflections,  and  Knudsen  layer. 

The  analysis  of  simulation  results  presents  the  characteristics  of 
the  laser  cutting  of  anode  for  lithium-ion  batteries.  Characteristics 
of  penetration  depth  and  absorptivity  change  significantly  at  the 
point  where  the  deep  penetration  hole  reaches  the  material  inter¬ 
face  due  mainly  to  the  composition  change  at  the  material  inter¬ 
face.  Furthermore,  temperature  distribution,  melt  pool  geometry, 
melt  pool  flow,  and  composition  distribution  are  examined. 

At  the  beginning  of  the  simulation,  the  L/V  interface  shows 
a  smooth  and  clean  surface  due  to  the  sublimation  characteristics 
of  graphite.  When  the  deep  penetration  hole  reaches  the  material 
interface,  a  melt  pool  starts  to  form,  develops  and  results  in  a  two- 
level  surface  due  mainly  to  the  sublimation  property  of  graphite 
and  the  melt  pool  flow  of  liquid  copper.  In  addition,  higher  values 
of  the  copper  concentration  are  observed  around  the  material 
interface.  When  the  deep  penetration  hole  is  becoming  deeper, 
the  copper  concentration  is  distributed  more  uniformly  and  the 
front  melt  pool  provides  strong  melt  pool  flow  in  the  copper 
region.  In  the  graphite  region,  the  copper  concentration  increases 
around  the  deformed  graphite  surface  and  melt  pool  flow  is 
observed  due  to  the  existence  of  the  liquid  copper.  Penetrations  on 
the  side  and  front  of  the  deep  penetration  hole  at  the  moment  of 
full  penetration  are  observed  since  uneven  melt  pool  shape  in  the 
copper  region  reflects  the  laser  beam  to  the  graphite  region.  The 
copper  concentration  is  relatively  low  around  the  bottom  of  the 
penetrated  hole.  Stronger  flow  is  observed  around  the  deformed 
graphite  surface. 

The  cut  surface,  kerf  width,  and  copper  composition  of  the 
experimentally  obtained  sample  are  measured  by  SEM  and  EDX 
analyses.  A  good  cut  surface  is  observed  and  kerf  widths  show 
a  discrepancy  with  simulation  results.  The  discrepancy  may  be 
caused  by  the  properties  of  the  graphite  sample  used  for  the 
experiments,  the  extended  copper  and  graphite  binary  diagram  in 
the  simulation,  or  the  lack  of  oxygen  reaction  and  air  suction 
included  in  the  simulation.  In  spite  of  these  discrepancies,  both 
experimentally  and  computationally  obtained  copper  concentra¬ 
tions  share  key  characteristics.  The  copper  concentration  starts  to 
increase  its  value  at  least  10  pm  above  the  material  interface  and 
increases  very  sharply  up  to  100%  copper  concentration  after  its 
initial  increase.  In  addition,  the  copper  concentration  decreases 
once  100%  copper  concentration  is  achieved. 

The  proposed  mathematical  model  of  laser  cutting  of  an  anode 
for  lithium-ion  batteries  can  be  utilized  to  predict  and  prevent  the 
defects  or  thermal  stress  observed  during  the  laser  cutting  process 
by  analyzing  the  melt  pool  and  temperature  distribution.  Further¬ 
more,  laser  parameters  can  be  optimized  through  this  model  to 
minimize  the  mass  transfer  of  copper  into  the  graphite  region. 
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Physical  constant 

R  gas  constant,  8.314  J  K^1  mol1 

a  Stefan-Boltzmann  constant,  5.67  x  10-8  W  m-2  K~4 
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incident  angle,  ° 

absorption  coefficient  of  copper 

absorption  coefficient  of  graphite 

liquid  constant-pressure  specific  heat,  J  kg-1  K-1 

solid  constant-pressure  specific  heat,  J  kg-1  K-1 

average  constant-pressure  mixture  specific  heat, 

J  kg"1  K"1 

average  constant-pressure  specific  heat  of  liquid, 

J  kg"1  K-1 

average  constant-pressure  specific  heat  of  solid, 

J  kg-1  K-1 

delta  function 

unit  vectors  in  xy,z  directions,  respectively 

emissivity 

liquid  mass  fraction 

solid  mass  fraction 

liquid  volume  fraction 

solid  volume  fraction 

force  function,  m  s-1 

speed  function  of  the  L/V  interface  due  to  evaporation, 
m  s_1 

enthalpy,  J  kg-1 
liquid  enthalpy,  J  kg-1 
solid  enthalpy,  J  kg-1 
smoothing  function 
incoming  ray 

thermal  conductivity,  W  m-1  K”1 
liquid  thermal  conductivity,  W  itT1  K_1 
solid  thermal  conductivity,  W  m_1  I<-1 
curvature 

latent  heat  of  fusion,  J  kg-1 
latent  heat  of  vaporization,  J  kg-1 
viscosity,  m2  s-1 
vapor  viscosity,  m2  s_1 
liquid  viscosity,  m2  s-1 

net  mass  loss  due  to  evaporation  per  area,  kg  s-1  m-2 

normal  vector 

pressure,  Pa 

density,  kg  m-3 

liquid  density,  kg  m~3 

solid  density,  kg  m~3 

the  actual  net  energy  influx  at  the  Liquid-Vapor  interface 
per  area,  J  s-1  m~2 

the  spatial  laser  beam  distribution  after  multiple 

reflections  per  area,  J  s_1  m~2 

energy  loss  due  to  evaporation  per  area,  J  s-1  m-2 

specific  gas  constant,  J  kg-1  K-1 

reflected  ray 

surface  tension,  N  itT1 

time,  s 

temperature,  I< 
temperature  at  surface,  K 
melting  temperature,  K 
boiling  temperature,  K 
temperature  at  vapor,  K 
fluid  velocity  vector,  m  s_1 
velocity  of  vapor,  m  s-1 
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x  spatial  vector 

(j)  level  set  function 


Abbreviations 

L/V  Liquid-Vapor  interface 

S/L  Solid-Liquid  interface 

HEV  Hybrid  Electric  Vehicle 

PHEV  Plug-in  Hybrid  Electric  Vehicle 

EV  Electric  Vehicle 
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